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Abstract 

We have performed a detailed Monte Carlo (MC) simulation for the Advanced Thin 
Ionization Calorimeter (ATIC) detector using the MC code FLUKA-2005 which is 
capable of simulating particles up to 10 PeV. The ATIC detector has completed two 
successful balloon flights from McMurdo, Antarctica lasting a total of more than 35 
days. ATIC is designed as a multiple, long duration balloon flight, investigation of 
the cosmic ray spectra from below 50 GeV to near 100 TeV total energy; using a fully 
active Bismuth Germanate(BGO) calorimeter. It is equipped with a large mosaic 
of silicon detector pixels capable of charge identification, and, for particle tracking, 
three projective layers of x-y scintillator hodoscopes, located above, in the middle 
and below a 0.75 nuclear interaction length graphite target. Our simulations are 
part of an analysis package of both nuclear (A) and energy dependences for different 
nuclei interacting in the ATIC detector. The MC simulates the response of different 
components of the detector such as the Si-matrix, the scintillator hodoscopes and 
the BGO calorimeter to various nuclei. We present comparisons of the FLUKA-2005 
MC calculations with GEANT calculations and with the ATIC CERN data. 
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1 Introduction 


High energy cosmic rays, particles of matter from distant regions of our galaxy 
and possibly elsewhere in the universe, can be directly observed by balloon- 
borne or space-based experiments. Understanding the cosmic ray composition 
and energy spectra can provide clues to the origin of cosmic rays, their acceler- 
ation mechanism, and their propagation through the galactic and intergalactic 
media. Cosmic rays contain all natural elements from proton to nuclei beyond 
nickel. The all-particle spectrum of cosmic rays (where species can not be 
distinguished) obeys a power law dependence of the observed flux on energy 
from 10 GeV up to the highest energies measured by ground-based air shower 
array experiments, over 10 2O eV. From these measurements, it is known that 
the spectrum is somewhat steeper above 10 16 eV than it is below 10 14 eV. This 
steepening is called the ’’knee”. It is not known whether this knee is due to: 

• a change in the propagation of cosmic rays; 

• a change in the acceleration mechanism of cosmic rays; 

• a change in the composition of the cosmic rays; 

The Advanced Thin Ionization Calorimeter (ATIC) is a balloon-borne experi- 
ment designed to investigate the elemental spectra of protons and other nuclei 
up to Z = 26, for energies from 50 GeV to 100 TeV, with a statistical accuracy 
of better than 30% for protons at the highest energy. (1) The ATIC balloon 
payload was launched to an altitude of 36 km above sea level from McMurdo, 
Antarctica in December 2000 and December 2002. These two flights lasted 
more than 35 days. 

In this paper, we discuss a data analysis tool we have developed based on 
FLUKA-2005. This code provides the capability of simulating high energy 
intercations up to 100 TeV, the highest energy measured by the ATIC detector. 


* Corresponding author. 

Email address: phguna@feynman.pfiys.subr.edu (R. M. Gunasingha). 

URL: http://suphys.phys.subr.edu (R. M. Gunasingha). 

1 Present address: Purple Mountain Observatory, Chinese Academy of Sciences, 
China 


2 




2 Scientific Objectives 


The ATIC detector focuses on direct measurements of high energy cosmic rays. 
The main scientific objectives of the ATIC experiment are: 

• to measure the H and He energy spectrum from ss 10 10 eV to ~ 10 14 eV 
with a single instrument; 

• to investigate the spectral differences between protons and helium; 

• to determine if the proton spectrum has a bend between 2 TeV and sw 40 
TeV; 

• to accurately measure the H/He ratio as a function of energy; 

• to determine if the spectra of heavy nuclei have the same indices; 

• to search for high-energy electrons and neutrals; 


3 ATIC Detector 


The ATIC detector shown in figure 1 is comprised of a fully active 50cm x 
50cm x 20cm deep BGO calorimeter, preceded by three 10 cm thick graphite 
targets sandwiched between panels of scintillator hodoscopes, having an out- 
wardly projective angle of 24 degrees. It is equipped at the top with a large 
mosaic of silicon detector pixels capable of charge identification. The BGO 
crystals of dimensions '2.5cm x 2.5 cm x 25.0cm long are arranged horizontally 
in 8 layers, with the long axis aligned alternately along the X and Y axes. 
The vertical direction is along the Z axis. There are 320 BGO crystals with 
40 crystals in each layer. The construction and commissioning of the detector, 
the data acquisition system and the software development for the detector, 
and detector calibration are discussed elsewhere (1). 



Fig. 1. ATIC detector, showing the Si-matrix, the three scintillator hodoscopes SI, 
S2, and S3, the graphite targets and the BGO calorimeter. 
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4 GEANT-3.21 ATIC Monte Carlo 


We have written a MC code for the ATIC detector using GEANT-3.21, that 
generates events with energies less than 20 TeV. Note that GEANT-3.21 is 
presently limited to energies below 20 TeV for hadrons. Figure 1 was generated 
by GEANT-3.21 using a geometry input file for the ATIC detector. 

In 1999, the ATIC detector was calibrated using test beams at CERN. The 
detector was illuminated with 375 GeV protons, 300 GeV electrons, 150 GeV 
electrons and 150 GeV pions. Figure 2 shows the energy spectra from data 
and MC simulations for 375 GeV protons. The data and the simulated results 
agreed well as expected. Figures 3 and 4 show the comparison of data for 150 
GeV protons and pions with our GEANT3.21 MC calculations. Note the MC 
calculations have been performed for the same number of events as the data. 



Fig. 2. Comparison of CERN 375 GeV proton data with GEANT. 



Fig. 3. Comparison of CERN 150 GeV proton data with GEANT. 

GEANT is widely used in high energy physics for the simulation of primary 
particles such as p, e and fx. When GEANT4 was introduced the support for 
GEANT-3.21 became limited. GEANT4 is evolving very rapidly but it still 
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Fig. 4. Comparison of CERN 150 GeV pion data with GEANT. 

has some limitations for simulations at very high energies ( TeV to PeV range) 
and cannot provide simulations of heavy ions. 

One of the primary goals of ATIC is to determine the indices for the energy 
spectrum of heavy ions from He, C, N, 0 to Fe up to very high energies. For 
this purpose we had to find a reliable simulation code to simulate heavy ions. 
After a thorough search, we decided to use FLUKA. 


5 Description of the FLUKA Code 


FLUKA which stands for “FLUktuierende KAskade” is a MC code able to 
simulate the transport and interaction of hadronic and electromagnetic par- 
ticles on any target material. It is a multi-particle, multipurpose code and is 
used in many different fields such as medical physics, astrophysics and high 
energy physics (2). It has been tested and verified with real data. For light 
nuclei, it uses the Fermi break model. Above 5 GeV/nucleon, the Dual Parton 
Model (DPM) is used. FLUKA can treat combinatorial geometry. It is widely 
used for studies related to both basic research and to applications in radiation 
propagation, dosimetry, cosmic ray and accelerator physics. Followings are a 
few areas to which modern FLUKA has been successfully applied (3). 

• Accelerators and shielding - CERN, SLAC; 

• Background and radiation damage studies- ATLAS and LHC experiments; 

• Dosimetry and radiotherapy- NASA; 

• Calorimetry - ATLAS test beams and ICARUS neutrino experiment; 

One of the main reasons for using FLUKA for ATIC is the ability to trans- 
port heavy ions at very high energies. FLUKA is designed to transport heavy 
ions and their interaction products from 5 GeV/nucleon to very high ener- 
gies, approximately 10 PeV, using a DPMJET-3 code which is a high energy 
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hadron-hadron, hadron-nucleus, nucleus-nucleus interaction model. It can run 
on various platforms such as Linux, Unix and HP. The validity of the models 
implemented in FLUKA have been benchmarked against several experimental 
data, from accelerators to cosmic ray showers (2). The FLUKA code is main- 
tained and updated regularly. The FLUKA discussion archive provides much 
support for the FLUKA user community (4). 


6 THe ATIC-FLUKA Code 


We implemented the ATIC geometry in the FLUKA code. Events were ran- 
domly generated in position, direction and energy above the Si-matrix plane. 
An energy dependence of E~ 2 ' 7 was assumed. Figure 5 shows a 150 GeV pro- 
ton vertically incident on ATIC very close to the center on top of the detector. 


150 GeV proton ATIC-FLUKA 


60 



X(cm) 

Fig. 5. A 150 GeV proton interacting in the ATIC detector 


We have compared the energy deposition in the ATIC BGO calorimeter using 
GEANT-3.21 and FLUKA. Figure 6 shows this comparison. Figure 7 shows 
the comparison of energy deposition by 1 TeV protons in ATIC. Figure 8 
shows the energy distributions obtained from FLUKA and CERN data for 
375 GeV protons incident on ATIC. Although at lower enrgies the agreemnet 
between GEANT and FLUKA and CERN data are good at higher energy, 
namely 1 TeV, there is a 10% descrepancy between GEANT and FLUKA. 
Because there is no data avaible at 1.0 TeV for ATIC, it is not clear which 
simulation is better. 
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I 1 50 GeV proton FlUKA and GEANT I 



Fig. 6. Comparison of the deposited energy spectra in the BGO calorimeter for 150 
GeV protons from GEANT and FLUKA calculations 



Fig. 7. Comparison of the deposited energy spectra in the BGO Calorimeter for 1.0 
TeV protons from GEANT and FLUKA calculations 



Fig. 8. Comparison of the deposited energy spectra in the BGO calorimeter for 375 
GeV protons from CERN data and FLUKA calculations 

7 Heavy Ion Simulation 


One of the primary goals of ATIC is to determine the indices for the energy 
spectra of nuclei such as He,C,N,0 up to Fe. For this purpose, we implemented 
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the DPMJET-3 interface to the FLUKA input card file. Figure 9 shows a 12 
TeV 12 C nucleus vertically incident on ATIC close to the center of the detector. 


1 TeV per nucleon Carbon-1 2 ATIC-FLUKA 



-60 -40 -20 0 20 40 60 


X(cm) 

Fig. 9. A 12 TeV 12 C nucleus interacting in the ATIC detector 

Figure 10 shows an event display for a 40 TeV 14 N nucleus incident on ATIC 
with a randomly chosen position and direction. Figure 11 shows a 64 TeV 16 0 
nucleus entering the top of the detector with a randomly chosen position and 
direction. 


Randomly generated 40 TeV Nitrogen-14 event ATiC-FLUKA 



-60 -40 -20 0 20 40 60 


X(cm) 

Fig. 10. A 40 TeV U N nucleus interacting in the ATIC detector 
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4 TaV par nucleon (Total E 64 TaV) randomly generated Oxygen-16 



X(cm) 

Fig. 11. A 64 TeV 16 O nucleus interacting in the ATIC detector 

8 Event Generation and Reconstruction 


An important goal of the simulation for the ATIC detector is to calculate 
efficiencies for various cuts applied to the data, so that an absolute cosmic ray 
flux can be obtained. The spectra of cosmic rays is well described by an inverse 
power law in energy with differential flux given by || = CE~ 2 70 . We have 
written a user code to generate events according to the above power law with 
energies ranging from 10 GeV to 100 TeV. The position is randomly chosen 
to be between -40 cm and +40 cm in both x and y. All events were generated 
in a plane at 10 cm above the top of the detector. The spectra of generated 
events is shown in figure 12. 



Fig. 12. 10 GeV to 100 TeV randomly generated protons on ATIC using FLUKA 
We ran the ATIC-FLIKA MC for 1,000,000 protons with energies ranging 
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from 50 GeV to 100 TeV. After applying a fiducial volume cut on the particle 
trajectory 2.5 cm from the edge of the calorimeter in x and y, 250,000 events 
remained. Figure 13 shows the distribution of the energy deposited in the BGO 
calorimeter for the events passing the fiducial volume cut. The resulting MC 
generated data are in the same format as the actual data and can be analyzed 
in the same manner. 



Fig. 13. Energy deposited in the BGO calorimeter by events staifying the fiducial 
volume cut 

To insure the integrity of trajectory reconstruction, we required that every 
BGO plane satisfy an energy threshold selection of 1.0 GeV. The energy dis- 
tribution of such events is shown in figure 14. 



Fig. 14. The energy of the events with a mininum deposited energy of 1.0 GeV in 
each layer of BGO. 


The reconstruction of the trajectory was obtained from a x 2 fit using the 
weighted pulse-heights in the BGO crystals. The trajectory was them improved 
using the position of the hit pixel in the Si-matrix plane as well as any hit 
scintillator along its intersection with the hodoscope. To check the goodness of 
the fit, we extended the trajectory back to the plane where the MC events were 
generated. The distance between the generated points and the reconstructed 
points is shown in Figure 15. 
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Fig. 15. The distance between the reconstructed positions and the generated posi- 
tions on a plane 10 cm above the Si-matrix 

For the n C simulation, a total of 550,000 events out of 1,380,000 events passed 
the above-mentioned fiducial volume cut. Figures 16, 17, and 18 show gener- 
ated, accepted and threshold selected energies, respectively. Figure 19 shows 
the distance between reconstructed and generated positions. 



io 3 
10 2 
10 
1 

10- 1 1 10 10 3 10 4 10 5 

Fig. 16. 50 GeV to 100 TeV randomly generated 12 C events 



9 Conclusion 


A simulation code was developed for the ATIC detector using FLUKA-2005. 
This code can be used to simulate ATIC events for very high energy primary 
particles (p,e, and /i etc.) and heavy ions. The FLUKA-2005 code provides 
simulation of particles to the maximum energy « 100 TeV measured by ATIC. 
We compared MC data from standard GEANT-3.21 and FLUKA below 20 
TeV. The agreement is reasonable. This ATIC-FLUKA MC code provides a 
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Fig. 17. Energy deposited in the BGO calorimeter by 12 C events staifying the fiducial 
volume cut 



nr 1 i 10 ^ 10 s 10 4 10 s 


Fig. 18. The energy deposition of 12 C events where each layer has more than 1.0 
GeV energy deposition. 


pivotal role for the analysis of the ATIG data and will be a part of the analysis 
package for heavy ions at high energies, hot achievable by the GEANT-3.21 


MC. 
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Fig. 19. The distance between the generated position and the reconstructed positions 

for 50 GeV to 100 TeV randomly generated 12 C events. 
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